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Abstract 

An SU(4) extrapolation of the chiral unitary theory in coupled channels done to study the scalar 
mesons in the charm sector is extended to produce results in finite volume. The theory in the 
infinite volume produces dynamically the -D s *o(2317) resonance by means of the coupled channels 
KD, nD s . Energy levels in the finite box are evaluated and, assuming that they would correspond 
to lattice results, the inverse problem of determining the bound states and phase shifts in the 
infinite volume from the lattice data is addressed. We observe that it is possible to obtain accurate 
KD phase shifts and the position of the -D s *o(2317) state, but it requires the explicit consideration 
of the two coupled channels in the analysis if one goes close to the r\D s threshold. We also show 
that the finite volume spectra look rather different in case the D s *q(2317) is a composite state of 
the two mesons, or if it corresponds to a non molecular state with a small overlap with the two 
meson system. We then show that a careful analysis of the finite volume data can shed some light 
on the nature of the -D s *o(2317) resonance as a KD molecule or otherwise. 
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I. INTRODUCTION 



One of the challenges of lattice QCD calculations is to determine spectra of mesons and 
baryons. Although problems stem when addressing excited states which have decay channels, 
steady progress is being done in the field [l T2] . In the case of resonances, the identification 
of the states obtained in the box with those in the infinite volume is not automatic since 
scattering states show up as discrete states in the box. The problem is formally solved when 
the states decay in one channel and the framework of Luscher allows to produce phase shifts 
in the infinite volume, and eventually investigate resonance properties, from the levels in 
the box [H, M ■ 



Recently, the Luscher ap pro ach has been generalized to the case of multi-channel scatter- 
ing. This was done in Refs. [l5(] on the basis of potential scattering theory, while the authors 
of Refs. (l6l . 17 ) use non-relativistic effective field theory (EFT) for this purpose. More re- 
cently the ideas of Ref . 17[ in the meson scalar sector and of Ref . [l8[ , where the case of the 



A-resonance in nN elastic scattering is addressed, have been followed up in Ref. [19f] , where 
a new and practical rederivation of the approach has been done, generalizing it to include 
fully relativistic effects. Furthermore, new methods of analysis of lattice data are presented 
for the case of two coupled channels and a thorough investigation of the propagation of errors 
is done. The results of Ref. [l9| were done for nir and KK scattering, addressing the issue 
of the / (600), / (980), and a (980) resonances. The "synthetic" lattice data were produced 
using the chiral unitary approach (where these resonances are dynamically generated) in 
a finite volume. A follow up of Ref. [l9] using input from the Juelich model [2(j, which 



requires discretization in momentum space, is done in Ref. [2lf] and further work concerning 
the k(800) scalar resonance is done in Ref. [22| . 

In the present case we face the problem of the scalar mesons in the charmed sector, 
concretely the strangeness 5 = 1, charm C = 1, where an extremely narrow resonance, 
the D s *o(2317), appears. This problem has been looked at in the literature studying the 
interaction of coupled meson-meson channels in this sector and the resonance is found to 



be dynamically generated, mostly from the interaction of the KD, r]D s channels |23H26 
An analysis of scalar form factors determined from lattice calculations (2?| finds a large 
KD scattering length hinting at the existence of a bound KD state which is also identified 
with the -D s *o(2317). Additional work is done along the lines of Refs. [23h26| in Ref. 28| 
including the NLO chiral Lagrangian for the interaction of pseudoscalar charmed mesons 
with Goldstone bosons j29|. The bound state of the -D<j*o(2317) is also found and the KD 
scattering length is evaluated. For other heavy meson-light meson channels, comparison of 
the scattering length with QCD lattice results is done and good agreement is found. 

In the present work we take the "phenomenological" potential of Ref. [26| and use the 
KD, r)D s channels to evaluate the energy levels in a box using periodic boundary conditions. 
In a second step, we take these "synthetic" lattice data and address the inverse problem of 
extracting from them the KD phase shifts in the infinite volume and the position of the 
D s * (2317), which appears as a bound state of the system. We show that with a few data 
from the lattice (about ten) corresponding to two levels of the box with various values of 
L (the box size), we can reproduce the results in the infinite volume, both for the energy 
of the bound state and the KD phase shifts, with a good accuracy. For this purpose we 
use the method with two channels developed in Ref. [19j, but generalize the problem to use 
two-meson loop functions in dimensional regularization rather than in the cut off method 
used in Ref. [l9|. We observe that the consideration of the two channels in the analysis is 
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necessary for an accurate reproduction of the data in the infinite volume. Further we show 
that using several instruments of analysis one can conclude from the synthetic data that the 
D s * (2317) state corresponds mostly to a bound state of KD and is not a genuine state in 
the sense used by Weinber g |30f] to decide that "the deuteron is not an elementary particle" 
or by the authors of Ref. [3l| to decide that the "a (980) and /o(980) are not elementary 
particles". This means that the analysis of future lattice data in this sector could bring 
an answer to this question in the case of the -D s * (2317) state. For this purpose we find 
out which data are needed and with which precision the energy levels must be evaluated. 
These results are very opportune since the Tokyo lattice group is addressing this problem 
at present 32| . 



II. FORMALISM 



A. The chiral unitary approach in infinite volume 



In the chiral unitary approach the scattering matrix in coupled channels is given by 



T=[1-VG]- 1 V (1) 

where V is the matrix for the transition potentials between the channels and G, a diagonal 
matrix, is in our case the loop function of two meson propagators, which is defined as 

r _. f d 4 q 1 1 ( ] 

% J (27r) 4 q 2 — m 2 + it (P — q) 2 — M 2 + it U 

where m and M are the masses of the two mesons and P the fourmomentum of the external 
meson-meson system. 

In dimensional regularization this integral is evaluated and gives for meson-meson systems 



33|,|3J 



Gi{s,mi,Mi, 



+ 



(An) 2 { 



. . , m 2 M 2 -m 2 + s^ M 2 

Oi(tl) + log -f + log — ^ 

fi z 2s rrq 

log - (Ml - m\) + 2ViQi(Vi)) + log (a + (M\ 



(3) 



+ 2^~sQi(^~s) 



log (-s + (Ml - m\) + 2 v ^g i ( v ^)) - log (-s - (Ml - m 2 ) + 2^0.^) 



where s = E 2 , with E the energy of the system in the center of mass frame, Qi the on shell 
momentum of the particles in the channel, \i a regularization scale and di(\L) a subtraction 
constant (note that there is only one independent parameter, because a change in [i can be 
absorbed into a$). 

In other works one uses regularization with a cut off in three momentum once the q° 
integration is analytically performed [35[ and one gets 
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(2tt)3 2^(0) £ 2 - Mg) + c^(g)) 2 + ie ' 



<*M = <J™ % j + P, uj'M = JM] + P (4) 



In Ref. 33] the equivalence of the two methods was established. 

We consider s-wave coupled-channels of KD and 7]D S with 1 = for the description of 
the _D s *o(2317) as a dynamically generated state. The third channel considered in Ref. [26], 
r] c D s , plays a negligible role and is not considered here. The pseudoscalar mesons K and D 
are represented as the doublet of isospin 

^ ' ^ (5) 




The isospin 1 = KD state is given by 

\KD,I = 0>= — -=\K + D° + K°D + > (6) 
V2 

The state of r)D s is already a state of I = 0. By naming the channels, 1 for KD and 2 
for i]D s , the transition potentials are given by [26J 

V 22 = -77tV^(- s + 2t - + 2m^ + 6m^ - 4m 2 ] 

^ 2 = fi /q 1 / , ^ U -t)-(3 + l)(s-u)-m 2 D - 3m 2 K + 2<] (?) 



where 



% \ mo* 



We study only the s-wave interaction and, hence, we must project t and u over s-wave, 
which is given by 

u = m\ + m\ — —(s + m\ — m\)(s + m\ — m\) (9) 

where rrii are the masses for the process 1, 2 — > 3, 4. The variable t projected over s-wave, t, 
can be obtained from the condition 

s + u + t = m x + m 2 + m 3 + m 4 (10) 

The use of Eq. (JTJ) with the two channels that we consider leads to a dynamically gen- 
erated state using the dimensionally regularized G D function with /i = 1500 MeV and the 
subtraction constant a = —1.26 at the energy of 2317 MeV, which we associate to the 
D s «o(2317) resonance. 
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B. The chiral unitary approach in a finite box 



When one wants to obtain the energy levels in the box, one replaces the G function by 
a G in Eq. ([Tj), where instead of integrating over the energy states of the infinite volume, 
with q — | q\ being a continuous variable, as in Eq. (JH), one sums over the discrete momenta 
allowed in a finite box of side L with periodic boundary conditions. We then have G = 
diag (Gi, G 2 ), where 



G; 



1 

Z 3 



l<?1<<?m; 

E 

q 



1 



— n, iifZ 

Lj 



(11) 



This is the procedure followed in Ref. [19| . 



Since the work of Ref. 26 



was done using G functions in dimensional regularization, we 
are going to adapt the finite volume formalism to the use of the dimensionally regularized 
G functions in the infinite volume. For this purpose let us write the G function of Eq. 



as 



G D (E) = b+ [ -^-I( q )+ lim / ^%-Uq)- lfm / ^\l a (q) (12) 



d 3 q 



d3( l r«W 



(2k) 



(2kY 



(2k) 3 ' 



where I(q) is the integrand of Eq. 



2uj (q) u< 3 (q) E 2 - ( Uj (q) + u' (q) ) 2 + it : 



(13) 



with Uj, u'j defined in Eq. (@J and I a (q) the asymptotic expression of I(q) when q goes to 
infinity. 



i 



2q 



1 



2q 2 -(2q) 2 + ie -Aq 3 



(14) 



The constant b in Eq. (TT2"j) is a remnant subtraction constant inherent to the regularization 
procedure. Thus we can write 



l<?1<M 



G D (E) = b+ I 7^3/(9)+ lim 

7 (27r) 3 9max->oo 



<? max 



^ J(g) + iln'^ 



(2k) 



8tt 2 



6 + lim 

<Jmax->°0 



d3q J(g) + iln /(/ - 



.<?<9n 



(2k 



(15) 



When we study the problem of the energy levels in the box we will substitute G D by G, 
replacing the integral by the discrete sum of Eq. ffTTl) up to the same q max , and since from 
Eq. f[T5"j) we have in the limit of q max — > oo 
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FIG. 1. Real part of the last two terms of the right hand side of Eq. (|17|) for the KD channel. 
The solid line indicates the average that we take between 1200 MeV and 2000 MeV for q max . The 
results correspond to a value of L = 2.4 m" 1 and E= 2317 MeV. 



6+ i in m= cD < £ >- / ^ m (16) 



<?<<?! 

we can write 



G(E) = G D (E) + lim 



"I <?max » 

27 ?'(*>- / 



9<<?n 



d 3 q 



(17) 



The three dimensional sum in Eq. (|T7|) can be reduced to one dimension considering the 



multiplicities of the cases having the same n [2ll. l36l| . The integral in Eq. i fPTT) has an 



37|] (see erratum). 



analytical form as shown in the appendix of Ref. 

When calculating the limit of q max going to infinity in Eq. (|T7j) we obtain oscillations 
which gradually vanish as q max goes to infinity. Yet, it is unnecessary to go to large values 
of q ma x, and performing an average for different q max values between 1200 MeV and 2000 
MeV one obtains a perfect convergence, as one can see in Fig. [TJ Note that the imaginary 
part of G D and of the integral in Eq. (1171) are identical and they cancel in the construction 
of G, which is a real function. 

The eigenenergies of the box correspond to energies that produce poles in the T matrix. 
Thus we search for these energies by looking for zeros of the determinant of 1 — VG 



det(l -VG) = 1- V n Gi - V 22 G 2 + {V n V, 



22 



V? 2 )G X G 2 = 



(18) 



In Fig. [2]we show the energy levels obtained for the box for different values of L. We show 
there the first 5 levels. We observe a smooth behavior of the levels as a function of L. The 
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FIG. 2. Energy levels as functions of the cubic box size L, derived from the coupled channels 
unitary approach of Ref. [26] and using G from Eq. (JTTJ) . 

lowest level converges for large L to the value of the energy of the bound state of the infinite 
volume case. This is a standard feature for the ground level bound state in the lattice, where 
the energy obtained is about the same as in the infinite volume for sufficiently large values 
of L. The first valuable information from our study is how large should L be to provide an 
accurate value of the energy. We find that with values of L = 3 m^ 1 there is already a good 
convergence to the value in the infinite volume. Yet, we would like to get more information 
from the lattice data, for instance scattering phase shifts of KD, and eventually i]D s , but 
we will restrict ourselves to just the first channel. We could also ask ourselves about the 
nature of the bound state found. Of course, in the present problem where the _D s *o(2317) 
appears dynamically generated we should get as an answer that it is indeed a bound state 
of KD with some admixture of rjD s , but if the lattice data were different than the levels 
obtained by us, the possibility exist that the answer would be different. We address these 
problems in section V. 



III. THE INVERSE PROBLEM OF GETTING PHASE SHIFTS FROM LATTICE 
DATA 



A. Results with two channels 



In this section we face the problem of getting bound states and phase shifts in the infinite 



volume from the energy levels obtained in the box using the two channel approach of Ref. [26 
which we would consider as "synthetic" lattice data. To accomplish this we need more 
information than just the lowest level, but we shall see that the first two levels shown in 
Fig. [2] already provide the necessary information to reproduce the problem in the infinite 
volume. 

In Ref. [3] several methods were suggested to solve the inverse problem, but we borrow 
here the one based on a fit to the data in terms of a potential suggested by the work of 
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Ref. H or Ref. [H-^. As we can see in Eq. ([7j), the potentials have a large constant part, 
some terms proportional to s and some terms inversely proportional to s. It is very easy 
to see that if one chooses a region of energies around a certain value of s, so, the inverse 
function of s can be expanded as a function of s — s to a good approximation. Choosing 
s = {m K + M D ) 2 then the ansatz of the following equation 



V, 



+ bij(s - (m K + M D ) 



(19) 



is a very accurate assumption. Comparing Eq. (fT9"j) with Eq. J7]), used in the chiral unitary 
approach of Ref. 26|| , we find 



a u = -155.101, b u = -3.732 x 10" 5 MeV~ 2 

a 12 = a 2i = -90.756, b 12 = b 21 = -3.361 x 10~ 5 MeV~ 2 (20) 
a 22 = -52.356, b 22 = 9.395 x 10" 7 MeV~ 2 



Certainly, this is not the unique option to face the inverse problem. In Ref. [19| it was 



already discussed that one could use a different parametrization and the consequences of 



it were discussed there. This problem has been studied in more detail in Ref. [22J where 
a systematic study of results obtained with different options of the potential is done. The 
conclusion of that work is that the results obtained with the simple form of the potential, 
Eq. (TT9~j) . are fairly good, but the freedom to use other potentials reverts into a somewhat 
larger uncertainty in the parameters of the resonances searched. 

We assume that the lattice studies provide us with ten eigenenergies corresponding to 
the first two levels of Fig. [2] for different values of L between 1.7 m~ x and 3.3 m~ l . We 
also assume that the levels are provided with an error of ±10 MeV, something achievable 
in present QCD lattice calculations. We make a best fit to the data assuming a potential 
as in Eq. (1191) . We look for the minimum \ 2 an d obtain a set of parameters for 0^,6^. 
Then we generate random sets of the parameters close to those of the minimum, such that 
y 2 is only increased below xLin + 1 ; a criteria that provides a band of fair statistical errors 
(l9l . 38| . With these values we generate the spectrum of Fig. E]by searching for the zeros 



of the determinant of 1 — VG. It might look like this is a tautology, since by using the 
same potential as in the original infinite volume problem and the same subtraction constant 
in G we should obtain for the parameters the same results as in the original potential (see 
Eq. ( 1201 ) and x m in = 0- And this is indeed the case. However, when the lattice data 
are provided to us we do not know which implicit regularization subtraction constant the 
lattice data would support. The inverse method only has a real value if the results that one 
obtains are independent of this subtraction constant. In one channel the improved Luscher 
approach of Ref. [l9| shows clearly that the results of physical quantities, like the phase 
shifts, are independent of the cut off chosen, for values of q max relatively larger than the 
on shell momenta of the particles (q m ax > 1-5 GeV in general), and has a defined limit for 



Qmax going to infinity. In the two coupled channels it was also found in Ref. [19] that the 
results were basically cut off independent, because a change in the re gula rization scheme 
should be absorbed by the fitted coefficients of the potential, Eq. ( TT91 . [39]. This is related 



to the basic feature of the renormalization group, with a trade off between the scale and the 
potential [4(| . This finding is what renders very valuable the approach of Ref. 19j to solve 



the inverse problem. The arbitrariness in the choice of the cut off in Ref. [19( translates here 
in the arbitrariness in the choice of the subtraction constant a in the analysis of the inverse 
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FIG. 3. Energy levels as functions of the cubic box size L, reconstructed from fits to the "data" 
of Fig. [2] using the potential of Eq. (|19p . The band corresponds to different choices of parameters 
within errors. 

problem, and we also find that the results are independent of this choice. This means that 
we choose an arbitrary subtraction constant in G in the analysis of the lattice data, and 
the same one in G D to reconstruct the phase shifts in the infinite volume. Like in the cut 
off method, we also take a range for the subtraction constant of a e [—1, —2], considered of 
natural size in Refs. 33l and 2(1 . 



As we mentioned above, in Fig. [3] we show now the results of the levels reconstructed from 
the best fits to the data, with a band corresponding to the random choices of parameters 
satisfying the condition that x 2 < + More interesting is to construct from these levels 
the KD phase shifts in the infinite volume. This can be seen in Fig. |H The phase shifts are 
evaluated from the Tn matrix element obtained from the Bethe-Salpeter equation, Eq. (pQ), 
using the potentials obtained in the fit and the continuum G D function. The normalization 



that we use is such that in one channel 35 



T{E) = , (21) 

pcot 8(p) — ip 

from where we determine the phase shifts in the infinite volume problem. 

In Fig. H]we can see the KD s-wave phase shifts for 1 = 0, 8$, reconstructed with our 
procedure for the different values of the parameters generated. As we can see, the agreement 
with the exact results is quite good, and we see how the errors in the determination of the 
lattice levels have propagated in the determination of the phase shifts. 



B. Results with one channel analysis 



In the present problem the 1]D S channel is far away from the KD threshold and the 
channel KD is more important than the r)D s for the bound state and for low energies of the 



KD scattering 26]. We may wonder whether a fit to the lattice data would be possible with 
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FIG. 4. Phase shifts for KD scattering derived from the coupled channels unitary approach of 
Ref. 26] (solid line). The band corresponds to the results obtained from the fits to the "data" of 



Fig. [5] using the potential of Eq. (fTUj) with two channels. 



only one channel. It is well known that the effects of far away channels which are not too 
relevant in a problem can be incorporated with modifications in the pot ential or G function 
of the main channel. This is sometimes done explicitly, like in Ref. |4l|, where an effective 
KN interaction is constructed that incorporates the effects of the 7rS channel, even if this 
latter channel is relevant. The reason for this is that if a second channel is far away from 
the region of energies studied, its effect in this region is fairly energy independent (changes 
in AE of order AE/(E2 — E), where E% is the threshold of this second channel and E the 
energy under investigation). Then the effect of this channel is easily accounted for by a 
slight modification in the subtraction constant of the first channel. Yet, as we approach the 
threshold of the second channel this trade off does not hold and the explicit need of that 
channel becomes manifest. 

For the one channel problem, as shown in Ref. the T-matrix in the infinite volume 



can be obtained for the energies which are eigenvalues of the box by 



T(E) = (V-\E) - G D (E)) = (G(E) - G D (E)) . (22) 

One can equally use the same procedure as done for two channels eliminating the Vn and 
V22 parts of the potential. The results with both methods are basically identical. By using 
the same subtraction constant as used to generate the spectra we obtain an = —191.954 
and 6n = —9.868 x 10 -5 MeV~ 2 for the best fit, similar to the value in the two channel 
case. Note, however, that as shown in Eq. f|2^|) the results for T are independent of this 
regularization constant. 

It is interesting to remark that given the structure of G in dimensional regularization, 
Eq. (|17p . in the difference of Eq. f[2"2"j) the function G D cancels identically and one finds 
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5. Phase shifts for KD scattering derived from the coupled channels unitary approach of 
26] (solid line). The band corresponds to using the fits to the "data" of Fig. [2] using the 



FIG. 
Ref. 

potential of Eq. ([TO]) with only the KD channel. 



T(E)- 1 



lim 

Qmax ^OO 



I Qmax » 

Li a , Jg<qn 



d 3 q 



(23) 



This result is the same as the one obtained in Ref. 3] starting with cut off regularization, 
and, as proved in Ref. 19| is nothing else than Liischer formula, except that Eq. f[2"31 keeps 
all the terms of the relativistic two body propagator. Indeed Lusher results have an implicit 
approximation since some terms of the real part of the two body propagator, which are 
exponentially suppressed in the physical region, are omitted (see section 2.2 and Appendix 
A of Ref. 0). 

In Fig. Owe can see the results for the phase shifts obtained with one channel. We observe 
that for low energies the results are very similar to those obtained with two channels. As 
the energy increases, the results with two channels approach better the exact results. 



IV. THE BOUND STATE AND ITS NATURE 



For the energy of the bound state we obtain also good results with the one and two 
channel methods. The numerical results are E = 2317 ± 5 MeV in both cases. 

It would be interesting to see if from the lattice data we could say something about the 
nature of the D s * (2317) state. Obviously with the "synthetic" lattice data which we have 
produced the answer is trivial since the state was dynamically generated in our approach. 
But if the real lattice results were different than those obtained here the answer might be 
different. 
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A. Dynamically generated states 



The problem has been solved for bound states close to threshold using the method of 
Weinberg, based on the knowledge of the scattering length and effective range 30j, and also 



31 



for resonances not far from a threshold 

that the deuteron was a composite state of neutron and proton 
the present problem has also been used in Ref. 



Using this method in Ref. [19| it was found 
The method, updated to 
where claims are made that knowing 
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the dependence of the KD scattering length on the K mass can help to learn about the 
nature of the -D<j*o(2317). 

We follow here a different approach by going to the root of the derivation of the results 



of Ref. [30|. These results stem from a sum rule that comes from the normalization to unity 
of the wave function of the bound state. A modern formulation of this sum rule can be seen 
in Ref. 43| and states that 



Ed 



dE 



-1 



(24) 



E=E a 



where gf are the residues of the Ta scattering matrices (coupling squared) at the pole of the 
bound state (a) and Ga are the propagator of the two particles of the corresponding channels 
(the loop function G D that we use here) . However, the normalization of the couplings and the 
G functions in Ref. 1431 correspond to a nonrelativistic Quantum mechanics formalism, while 
here we use a different normalization in the field theoretical approach (the correspondence 
can be found in Ref. (43[). It is easy to generalize this result for the field theoretical case as 



Y.9 



ds 



(25) 



E=E n 



where now the residues are defined as 



9i 



lim (s 



sr)Tu 



(26) 



with sr the energy squared of the bound state. A different derivation of this equation is 
done in Ref. 44| using the Ward Takahashi identity. For the case of one channel this is easy 
to see since 



T 



1 



V- 1 - G 

and if V is energy independent it follows immediately that 



1 



9 



dG/ds 



(27) 



(28) 



One can see from the derivation in Ref. [43J that each term in Eq. (1251) accounts (with 
reversed sign) for the probability of the bound state to be a bound state of the pair of 
particles of the channel considered. In the case that we had the coupling of the bound state 
to another hypothetical elementary particle outside the space of pair of particles considered, 
there would be an extra term in the sum of Eq. (TSoT) . -Z, where Z would account for the 
overlap of the bound state with this hypothetical genuine particle. The diversion of the sum 
of Eq. (121)1) with respect to -1 would indicate the amount of the bound state which cannot 
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be considered a bound state of the two particles. Actually, for the compositeness of a bound 
state in a single channel, the relation 

-,*f = l-Z (29) 

was derived with a field theoretical argument j45|. Here Z is the field renormalization 
constant for the genuine state introduced to the theory by hand. 

In our formulation, there is a caveat, since in the derivation we have assumed that V is 
energy independent, while we are now producing potentials with a moderate energy depen- 
dence. We have checked that if we use a potential independent of the energy with the two 
channels, the theorem of Eq. (1251) holds exactly in our calculation. Yet, with the energy 
dependence of Eq. ( 1T9"]) the sum rule provides a value of -0.8. We should not look at this 
as an indication that we have the coupling to a genuine state. We should first see what 
comes out from the consideration of the energy dependence which we know explicitly in our 
problem. Taking into account the e nerg y dependence of Eq. (I19p we find that approximately 
one should have now (see also Ref. [441 ]) 



i % \ [an + hi(s - (m D + m K ) 2 )] 2 ds 

The former equation is fulfilled up to the level of 3%, hence its application could tell us, 
if the diversion from -1 of the sum rule is appreciable, that there would be a substantial 
coupling to a genuine state in the present work. Due to the construction of the state in 
the present problem, this is obviously not the case, as one can see by the proximity to -1 of 
the sum in Eq. (!3Up . Furthermore, we observe that about 90 % of the sum rule comes from 
the KD state, indicating that we have largely a bound KD channel in our approach. The 
successful analysis of one channel would obviously indicate that the state is 100% a bound 
state of KD. This shows the difficulty that we might have quantizing with high precision 
the amount of KD in the bound state of -D s *o(2317) from the analysis of the lattice data, 
yet, the use of an extra level of the box would allow us to be more precise. 




B. The case of a genuine state 

On the other hand, there is the possibility that the data might be such that one could find 
a gross deviation from having KD as the large component of the _D s *o(2317) wave function, 
reflecting the fact that the D s * (2317) would be a genuine state, not generated by the KD 
interaction. For this purpose we have made a test introducing a different potential where a 
CDD pole (Castillejo, Dalitz, Dyson) [46] is introduced by hand. The potential in just one 
channel would now be 



V = V M + 9 ° dd (31) 

S ~ S C DD 

where Vm is assumed to be energy independent and gc DD , scdd ar e the parameters of the 
CDD pole. 
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FIG. 6. Energy levels in the box constructed from the potential of Eq. (|31|) that contains a CDD 
pole. 

Using Eqs. (l26l) and (l2?l) and taking into account that for the bound state we have 
V" 1 = G, we find now (see also Ref. (4o| ) 



2 dG = 1 7 

9 ds dv- 1 /dGy 1 

ds \ ds J 
1 



9cddG 



(32) 



, 2 dG 



+ 1 



dG 

which indeed shows that Z is a quantity between and 1 since — — < below the KD 

ds 

threshold. 

We have chosen a potential of the type of Eq. (l3~Ii) with Vm of the order of ten times 
smaller that the potential used for Vn, ^scdd corresponding to a 20 MeV above the mass of 
the _D s *o(2317) and qcdd = 3787.34 MeV such as to reproduce the bound state at the mass 
of the -D s *o(2317). We find that at the pole of this state, Eq. f[3"2"j) provides Z ~ 0.9, thus 
showing that the introduction of a CDD pole as in Eq. (I3T1) is a good method to account for 
a genuine state beyond the space of the two interacting particles. 

The levels in the box with this potential are shown in Fig. As we can see, the levels 
show quite a different behaviour than those found with the potential of the coupled channel 
unitary approach [26| shown in Fig. [2J It is clear that the determination of the levels with 
lattice calculations can easily differentiate between these two scenarios. 

In Fig. [7| we show the results with a fit to the first two levels of Fig. O obtained with the 
CDD pole structure for the potential. We have taken also ten points over the curves and 
assumed ±10 MeV errors. With the values of the parameters that fulfill x 2 < Xmin + 1 we 
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FIG. 7. Fit to the first two energy levels in the box of Fig. [H] constructed from the potential of 
Eq. (f3Tj) that contains a CDD pole. 



show in Fig. [H] the phase shifts that one would obtain. We observe that the phase shifts are 
quite different from those obtained in the two channels with the standard potential. 

At this point we would like to comment that the phase shifts obtained with the standard 
potential, Eq. (Il9l) . or with the CDD potential, Eq. are negative close to the KD 

threshold. According to Ref. [ItJ (see also Ref. (48|) this is related to the presence of the 
bound state and Levinson's theorem which states that the difference of the phase shift 
between zero energy and infinity, for one channel with ordinary potentials, is me, with n the 
number of bound states. This induces a tendency for the phase shifts to start decreasing 
from threshold. This tendency has been used in lattice calculations [47| with the help of 
Luscher's theorem to induce that one has indeed bound states in the infinite volume case. 
Note, however, that as found here, this feature alone would not tell us wheather the bound 
state is dynamically generated from KD or a genuine state. Furthermore, we have checked 
explicitly that a potential of the CDD type does not fulfill the Levinson theorem., since we 
find that the phase shift decreases moderately as we increase the energy from threshold but 
then it starts increasing and hence does not decrease to 5(0) — n (see Fig. [S]) 



C. Case of a composite state analyzed in terms of a CDD potential 

It is also interesting to see if we could get a fit to the levels in the box obtained with the 
coupled channels approach of Ref. 26[ , shown in Fig. [2J by using the potential of Eq. (IBTj) 



that contains a CDD pole. Considering these energy levels as "lattice data" we fit them 
with the potential of Eq. (1311) and find the results shown in Fig. |HJ As we can see, the fit 
to the data is quite good. It is also interesting to plot the KD phase shifts obtained with 
this potential. In Fig. [TU] we show the results for the phase shifts obtained with the CDD 
potential. We can see that the phase shifts obtained up to 2400 MeV coincide with those 
obtained with the two channel potential. After this, there is a large dispersion of the results 
and they divert significantly from those obtained with the two channel approach. Certainly, 
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FIG. 8. Phase shifts for KD from the fits of Fig. [7J using the potential of Eq. that contains a 
CDD pole. The solid line corresponds to the results obtained with the CDD potential of reference, 
and the band corresponds to the fits to the first two levels with a potential of the same type 
considered in Fig [7J 

the use of an additional "lattice" level would put big constraints in the phase shifts in that 
region helping us decide between the two options. 

However, it is interesting to analyse the results obtained: we find the CDD pole at the 
energy of about 2500 MeV, far away from the bound state of the _D s *o(2317) state. If one 
restricts oneself to low energies, the fit with the CDD would be acceptable. Yet, the fact that 
the CDD pole has appeared so far away from the energy of the -D s * (2317) state is telling us 
that the data do not want a CDD pole being responsible for this state. The CDD pole far 
away in this case simply generates a smooth energy dependent potential in the region of the 
low energies. The interesting thing is that if we calculate now Z from Eq. (1321) we find that 
Z ~ 0.15, indicating that the bound state is basically a KD bound state, with the precision 
that the limited low energy data provide. A more precise determination would require the 
consideration of an additional level of the box. 

D. Case of a genuine state analyzed with an ordinary potential 

It is also interesting to perform another test. Let us assume we had the lattice data which 
correspond to the levels generated with the CDD potential in Fig. [H]and we would like to fit 
them with the smooth potential of Eq. (TTTJ]) . The fit is bad, the \ 2 is now °f the order of 6 
and the fit to the data can be seen in Fig. [TTJ The quality of the fit is worse than in all the 
other fits that we have produced. Certainly, a better precision in the lattice data, and/or 
the addition of extra data at smaller box sizes would tell us that the fit is actually very bad. 
The addition of a third level would certainly help to disregard this solution. The exercise 
has also served to show the importance of the scattering data to determine the nature of a 
bound state. 

To finalize the discussion, we present in Fig. [T5] the results of the phase shifts that this 
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FIG. 10. Phase shifts from the fits with a CDD pole to the first two levels obtained with the 
coupled channels approach of Ref. [26j] shown in Fig. [9l The solid line corresponds to the results 
with the two channel analysis of Fig. SJ 



fit would generate. As we can see, the results are very different from those generated from 
the CDD potential, such that if more precision is demanded to the lattice data such that 
they can produce the phase shifts with accuracy, this type of fit would be easily ruled out. 
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FIG. 11. Fit to the levels produced with a CDD potential using the potential of Eq. (|19|) with one 
channel. 
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FIG. 12. Phase shifts derived from the fit with the potential of Eq. (|19p with one channel to the 
levels produced with the CDD potential (lower curve). The upper line is the phase shift obtained 
with the CDD potential shown for reference. 



V. CONCLUSIONS 



In summary we have addressed three problems in the present work. The first one is to 
use the unitary coupled channel method in finite volume, generating the levels in a finite 
size box as a function of L. The second problem is the inverse problem: we assume that the 
results obtained before would correspond to results given by lattice calculations. From them 
we would like to obtain bound states and scattering states in the case of the infinite volume. 
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We found that the two channel method proposed by us works well and is independent of the 
regularization scale used, which obviously we do not know when the lattice data are provided 
to us. If we stick to low energies of the KD system we also see that the one channel analysis 
is quite good, but as we approach the threshold of t]D s , the two channel analysis provides 
better phase shifts. The third issue faced in the work is whether we could say something 
from the lattice results concerning the nature of the bound state obtained. We found that 
indeed it was possible to induce some information on the nature of this state as basically a 
molecule of KD, by inspecting the spectra in finite volume. However we needed two levels 
to get an answer to this problem. We carried out a useful exercise looking for the different 
alternatives in the analysis. We showed that if the bound state is a composite state of 
KD, the analysis using an ordinary potential and the Weinberg condition provided Z ~ 
indicating the compositeness of the state. Alternatively we tried to analyze the lattice data 
in terms of a CDD potential that could accommodate a non molecular state and we could 
find a good fit to the data but with an artificial pole at an energy very far away, such that 
the potential was essentially constant in the region of interest and the Z value for overlap 
with a genuine state was small, compatible with zero within the precision that the data 
allow. 

Similarly we played the exercise of assuming that the lattice data would correspond to a 
CDD potential that represents mostly a genuine state, with Z ~ 0.9. The analysis with a 
CDD potential provided a good fit to the data and returned this value of Z. On the contrary, 
attempts to fit these data with an ordinary potential failed, leading to unacceptable values 
of the x 2 - 

All these exercises show that the data provided by QCD lattice calculations on energy 
levels in a box contain the information to decide on the nature of the bound state D s *o(2317). 
With the amount of data and precision suggested here one can obtain Z, the probability to 
have a genuine state, with a precision of about of about AZ ~ 0.1. This is already quite 
a good precision at a time where the discussion goes on around the nature of some states 
at a qualitative level. Certainly more information and better precision on the lattice results 
could improve this error if desired. 
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